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1.  Introducti on 

The  class  of  bivariate  distributions  in  which  we  are  interested 

in  this  paper  are  those  which  are  positive  quadrant  dependent.  To  be  more 

specific,  let  8  be  the  Borel  a-field  on  the  real  line,  R  and  8x8  be 

2 

the  Borel  o-field  on  the  two-dimensional  euclidean  space,  R  .  If  y  is 
a  probability  measure  on  BxS  ,  we  denote  the  first  and  second  marginals 

of  U  by  and  U2  respectively.  The  probability  measure  y  is 

said  to  be  positive  quadrant  dependent  if 

y{tc,«)  X  |d,*)}  _>  yj^{{c,*)  }y2{id,'») } 

for  every  c,d  in  R.  The  above  notion  can  be  restated  as  follows.  Let 
X  and  Y  be  two  real  random  variables  with  some  joint  distribution  function 
F.  X  and  Y  are  said  to  be  positive  quadrant  dependent  if  \ 

“3  ■ 

Pr(X  ^  X,  Y  ^  y)  PrCX  >  x)  Pr(Y  >  y)  Q 

□ 

for  all  x,y  in  R.  See  Lehmemn  (1966)  for  this  and  related  notions  of _ 


dependence . 


By . 

Diit  tbulion/ 
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We  want  to  look  at  the  above  notion  of  dependence  from  a  global 
point  of  view.  Choquet's  theorem  is  in  the  background  of  the  ensuing  discussion. 
Suppose  M  is  a  compact  convex  subset  of  a  locally  convex  linear  topological 
space.  Then  by  Krein-Milman  theorem,  there  is  atleast  one  extreme  point  of  M. 
Choquet's  theorem  represents  every  member  of  M  as  a  mixture  (in  some  sense) 
of  the  extreme  points  of  M.  For  a  full  discussion  on  Choquet's  theorem, 
see  Phelps  (1966) .  We  want  to  piursue  the  same  line  of  tack  in  studying  the 
notion  of  positive  quadrant  dependence.  Accordingly,  let  denote  ^ 

the  collection  of^bivariate  positive  quadrant  dependent  distribtuions. 

We  equip  with  the  topology  of  weak  convergence.  Suppose  (i) 

is  compact  and  (ii)  is  convex.  Then  in  accordance  with  Choquet's  ^ 

theorem  we  can  write  every  member  of  *  mixture  of  extreme  points 

of  There  are  certain  common  properties  of  extreme  point  distributions 

which  are  preserved  under  mixtures  and  a  good  understanding  of  the  set  of 

all  extreme  point  distributions  would  provide  an  insight  into  the  structure 

of  a  general  bivariate  positive  quadrant  dependent  distribution. 

Ie 

But,  the  set  is  neither  compact  nor  convex.  Non-coi»actness  ^ 

of  the  set  is  not  surprising.  The  malaise  stems  from  the  non~compactness 

of  the  class  of  all  bivariate  distributions.  The  following  bivariate 

r 

distributions  are  positive  quadremt  dependent. 
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1  1/4 

1/4 

1 

2  1/4 

1/4 

_ 1 

1 

1 

2 

1 

2/9 

1/9 

2 

( 

4/9 

2/9 

1 

I  1 

But  (l/2)y  +  (1/2) X  is  not  positive  quadrant  dependent. 


•***.•  •  k** 
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In  the  next  section,  we  identify  some  natural  subsets  of  W 
which  are  compact  and  convex.  In  Section  3,  we  discuss  a  method  of  enumerating  I-  . 

I 

W-- 

the  extreme  points  of  these  compact  convex  sets .  In  Section  4 ,  an  application  >  v, 
of  the  results  of  Sections  3  and  4  is  discussed  in  detail  in  the  context 
of  contingency  tables. 


2.  Comoact  convex  subsets  of  M _ . 

- - - ?QD 


Let  X  and  v  be  two  fixed  probability  measures  on  the  Borel 
Q-field  B  of  the  real  line.  Let  .M(X,v)  be  the  set  of  all  bivariate  *-V.. 

distributions  y  with  the  first  and  second  marginals  X  and  v  ;1-. 

•  *1*  a 

respectively,  i.e., 

.M{X,v)  =  {  y;  y  is  a  probability  measure  on  BxB  ,  y  »  X  .v.v 

and  y.  *  v}  .  .'•.'■i 

Let  the  set  of  all  bivariate  positive  quadrant  dependent 

distributions  with  the  first  and  second  marginals  X  and  v  respectively,  i.e 

Mpgj^(X,v)  «  {y  ;  u  is  positive  quadrant  dependent,  y^^  »  X 

and  y^  -  v}. 

These  sets  are  topologically  and  functional  analvticallv  are  nice  •  j 
as  the  following  result  exemplifies. 

Theorem  2.1  M(X,v)  and  M _ (X,v)  are  comoact  convex  sets.  . 

-  PQD 

Proof.  See  BhasKara  Rao.  Krishnaiah  and  Subramanvam  (1985). 


Thus  we  obtain  a  decomoosition  of  into  conroact  convex  subsets  v'v'i 

PQD 


as  follows. 


u  u  ■v'’''’ 

A  V 


r' 


3.  Extreae  points 

Extreme  points  of  the  cossact  convex  set  have  been 

characterized  by  Lindenstrauss  (1965)  and  Douglas  (1964) .  The  following 
is  their  result. 

Theorem  3.1  A  probability  measure  U  in  M(X,v)  is  an  extreme  point  of 

M(X,y)  if  and  only  if  the  collection  of  all  functions  h  from  RxR  to  R 

which  are  of  the  form  h(x,y)  «  f (x)  +  g(y)  ,  x,y  in  R  for  some  f  and 

integrable  with  respect  to  X  and  v  respectively 
g  real -valued  integrable  functions  on  the  real  line^i^is  dense  in 

(RxR,  8x8,  u) . 

Proof .  See  Douglas  (1964,  Theorem  1,  p.243).  See  also  Lindenstrauss  (1965, 
Theorem  1,  p.379). 

we  generalize  the  above  result.  We  operate  in  the  more  general 
realm  of  probability  triplets.  Let  (X,C,X)  and  (Y,V,v)  be  two 
probability  spaces.  Let  M(X,v)  be  the  collection  of  all  probability 
measures  u  on  the  product  a-field  CxV  such  that  the  first  marginal  of 
y  is  X  and  the  second  marginal  of  y  is  v.  This  notation  agrees 
with  the  notation  introduced  above  in  the  context  of  probability  measures 
on  the  real  line.  Let 


where  I,  denotes  the  indicator  function  of  the  set  A.  Let  F  be  the 
A 


linear  manifold  spanned  by  F  . 


'Pieoroni  3.2  A  probability  measure  U  on  CxP  is  an  extreme  point  of 
if  and  only  if  F  is  dense  in  L^^CXxY,  CxP,u)- 

Proof.  The  proof  here  is  modelled  on  the  proof  of  Theorem  1  of  Douglas  (1964 
p.243).  Suppose  U  is  not  an  extreme  point  of  M(X,v).  Let  y  *  +  n)/2 

for  some  ?  and  h  in  M(X,V).  This  implies  that  2y  ^  h  ^  Of  and  by 
Radon-Nikodym  theorem,  there  exists  a  function  h  in  L^CXxY,  Cxt7,y)  such 
that 

^  •  h  a.e.  [y]  ,  and 

1  -  h  0  a.e.  [y]  . 

The  function  1  -  h  is  orthogonal  to  F  ,  i . e . , 

/f(l  -  h)  dy  -  /f  dy  -  /fh  dy 

-  /f  dy  -  /f  dn  »  /f  dy  -  ff  dy 

*  O  for  every  f  in  F. 

This  clearly  demonstrates  that  F  is  not  dense  in  Lj^(XxY,  CxP,y)  .  For, 
if  F  were  to  be  dense,  since  the  dual  L*  (XxY,  CxD,]i)  of  L^^CXxY,  CxC,y) 

is  L  (XxY,  CxV,\i)  ,  the  linear  functional  induced  by  1  -  h  on  L,  vanishes 
identically  on  F  would  imply  that  1  -  h  »  O  a.e.Qi]  .This  contradiction 
proves  the  assertion  made  above. 


Conversely,  suppose  F  is  not  dense  in  L^(XxY,  C;d),U).  Then 

there  exists  an  essentially  non-zero  function  h  in  L  (XxY,  Cid^/U) 

00 

orthogonal  to  F.  Set 

n  (E)  •  /  h  du  for  E  in  Cad?  . 

E 

Since  /  hf  «  O  for  every  f  in  F,  we  have,  in  particulatr,  /  h  du 
This  implies  that  r.  (XxY)  »  O-  Define 


^(E)  -  /  (1  +  h/llhll^)  du  ,  and 

E 

U(E)  -  /  (1  -  h/l|h|l^)  du  for  E  in  CxP  , 

E 

where  [[^Mao  L^-norm  of  h.  Since  1  +  (h/jlhjl^)  ^  O  a.e.  |V|  * 

the  set  functions  ^  and  u  are  non-negative.  In  fact,  s,  u  £  M(a,v). 

We  can  now  write  u  *  *  u)/2.  This  completes  the  proof. 


We  now  exhibit  some  extreme  points  of  M(a,v}  using  measure 
preserving  transformations.  Let  T  be  a  measure  preserving  transformation 
from  X  to  Y  preserving  the  measures  X  and  v,  i.e.,  T  ^(D)  z  C  and 
X(T  ^(D))  =  v(D)  for  every.  D  in  V.  Let  G  be  the  graph  of  T,  i.e., 

G  ■  { (x,  Tx)  ;  X  e  X}. 

Let  be  the  projection  map  from  XxY  onto  X,  i.e.,  Pj^(x,y)  »  x  for 

every  {x,y)  in  XxY.  We  claim  that  for  every  E  in  x  ,  P^ (E  H  g)  is 
available  in  C.For,  let  E  *  {s  e  CxV  ;  P^ (£  O  G)£  C}.  One  can  check  that 


the  rectangle  sets  in  are  availeible  in  E ,  E  is  closed  under 

conplenentation  and  countable  unions.  Hence  E  *  CxO.  Define  a  set 
function  y  on  as  follows. 

y  (E)  ■  X  (E  n  G) )  for  E  in  . 

we  now  check  that  y  is  a  probability  measure  on  with  maurginals 

X  and  V. 

y (Cxy)  -  X (P^ ( (CxY) O  G) ) 

=  X{cn  t”^Y)  =  X(C)  for  C  in  C. 

y(XxD)  =  X(?j^{(XxD)n  G))  for  D  in  V. 

»  X(XfXT"^D)  =  X(t'^D) 

*  v(D)  ,  since  T  is  measure  preserving. 

The  assertion  that  y  is  a  probability  measure  is  obvious. 

We  now  claim  that  G  is  a  thick  subset  of  XxY  under  y,  i.e. , 
the  outer  measure  of  G,  y* (G)  »  1.  For,  if  S  is  any  set  in  CxV 
containing  G,  then  P^^CEO  G)  »  P^(G)  *  X.  Loosely  speaking,  the 
measure  y  is  concentrated  on  the  graph  of  T. 

Finally,  we  assert  that  y  is  an  extreme  point  of  ,M(X,v).  Suppose 
y  *  {5  +  u)/2  for  some  5  and  u  in  M(X,v).  It  is  obvious  that 

5*(G)  ■  1  “  y*  (G) .  Fxirther,  for  any  E  in  CxV,  5*(E  H  G)  =  ;(E) 

and  ’j*{EftG)  '»  U(E).  If  CxD  e  CxV,  (CxD)  O  G  C(C  ^  T  ^D)  x  Y.  Consequently, 


C(CxD)  <  ;((CAT  D)  xY)  »  X(CnT  D)  -  p(CxD), 

and 

U(CxD)  £  \J(CnT"S)  X  Y)  »  A(cnT”^D)  =  U(CxD). 

Hence  ij(CxD)  =  ^(CxD)  *  u(CxD)  for  every  C  in  C  and  D  in  V. 

Therefore,  v:  *  ^  *  u.  This  completes  the  proof. 

We  now  want  to  discuss  characterizing  the  extreme  points  of 

where  X  v  fixed  probability  measures  on  the  Borel 

(j- field  S  of  the  real  line.  One  can  prove  that  the  product  measure  y  x  v 

is  always  an  extreme  point  of  •  See  Subramanyam  and  Bhaskara 

Rao  (1983).  But  this  is  in  total  contrast  with  the  case  of  M(X,v)*  Th® 

product  measure  y  x  y  is  an  extreme  point  of  M(X»v)  if  and  only  if 

either  X  is  0-1  valued  or  y  is  0-1  valued.  See  Kemp  (1968,  p.l356). 

We  do  not  have  any  characterization  of  the  extreme  points  of  MpQpUrv) 

in  the  general  case.  Note  that  Hp^p(X»v)  is  a  compact  convex  subset 

of  SUXrj)-  It  is  easy  to  see  that  if  y  is  an  extreme  point  of  M(X,'5>) 

and  belongs  to  then  y  is  an  extreme  ooint  of  ■ 

**  PyD  ^ 

If  y  is  an  extreme  point  of  ,  y  need  not  be  an  extreme  point 

of  M(X»v5'  The  numbers  also  do  not  exhibit  a  pattern.  If  X  and  y 

s 

have  finite  support,  there  are  case  when  the  number  of  extreme  points  of 
iWp^p(X,v)  is  strictly  less  than  the  number  of  extreme  points  of  M(X,u), 
and  cases  the  opposite  happens. 


When  the  probability  measures  X  and  v  have  finite  support, 

the  subject  of  the  determination  of  the  extreme  points  of  M(X,v)  and 

of  M  (X,V)  have  a  distinct  flavour.  We  now  proceed  to  describe  a 
PQD 

method  of  determining  the  extreme  points  of  these  compact  convex  sets. 
Without  loss  of  generality,  assume  that  the  support  of  X  is  {l,2,...,m} 
and  that  of  v  is  {l,2,...,n}.  Let  p^  =  X{{i}),  i  *  l,2,...,in  and 
Q  =  V({j}),  j  =  l,2,...,n.  From  now  on,  we  use  the  suggestive  notation  : 

'j 


M(X,u)  »  iM (p^ ,p2 ' *  *  •  ^1 ' ^2 '  * ' "  ' 


and 


“  ‘^^QD<-°1''=2 . -V°-l'°-2 . V' 


Extreme  points  of  ^'^{pj^,P2' •  •  •  ‘ 


Any  bivariate  distribution  in  rP2»  •  •  • ' '^n^ 

can  be  described  in  the  form  of  a  matrix  ?  =  order  m  x  n 

with  the  following  properties. 


(i) 

>  0 

for 

all 

i 

and  j . 

(ii) 

m 

Z 

i-1 

for 

j 

*  1,2,. 

(iii) 

n 

r 

j-1 

^3 

'  ?i 

for 

i 

-  1,2,. 

Definition  3.i  Let  P  =  (p^^)  be  any  matrix  in  ^(Pj^rPj/- •  •  •  * 

An  ordered  sequence  ^  ^  '  —  '  P,  i  J  **  said  to  be  a 

^1^1  2^2 

loop  in  P  if 


(i)  k  is  even. 


(ii)  i,  »  i., 


i-.  -  ^A> 


i.  .  *  i.,/ 

V—  I  1C 


(lii)  ^2  “  J3'  ^4  =  35' 

(iv)  the  pairs  ^  *  1/2,..., 3c  are  distinct. 

A  loop  in  P  is  said  to  be  positive  if  every  meusber  of  the  loop  is 
positive.  Support  of  P  is  the  collection  of  all  pairs  (i,j)  for  which 
p^^  is  positive. 


The  following  lemmas  are  helpful  in  characterizing  the  extreme 


points  of  M(o,  ,p.,,  —  ,0  ;q,  ,q.,»  —  ,q  )• 
12  ID  1  2  n 


Lemma  3.1  Every  row  (column)  of  a  matrix  P  in  M(p,  ,p-,  —  ,p  ;q,  ,q.,i 

'  12  Q  1  2 

contains  an  even  number  of  elements  of  any  loop  in  P. 


Proof.  Obvious. 


p^  ^  ,  p p.  ^  ,  p.  j  ;  ... 

^2^2  2^3  r-'r  r^r+1 


Lemma  3.2  If  p.  •  »  P.  •  an  infinite  sequence  of  elements  from  a 

‘^1^1  ^1^2^ 


a  given  matrix  P  with  the  property  that  any  two  consecutive  suffixes  are 
distinct,  then  we  can  find  a  loop  in  P  consisting  of  elements  from  the 
given  sequence. 


Proof .  Let  b  be  the  first  element  in  the  sequence  whose  suffix  agrees 
with  the  suffix  of  one  of  the  elements  a  preceding  b. 


Case  (i).  a  -  p.  .  and  b  =  p  .  .  Then  [p.  a  •  P-  u  i 

r^r  ^k^k+1 


p.  ^  ,  p.  .  J  is  a  loop  in  ?. 

1__  n  ^  e  1_ 


r-l^r-1  r-l^r 


Case  (ii) .  a  *  p  .  and  b  »  p.  .  .  Then  fp. 

^k  k+1  ^r  r+1  "  ^ 


j  '  *  i  j  ' 
k+l-'k+l  '^k+l-'k+2 


Pi  .  ,  Pi  .  J  is  a  loop  in  P. 
^r^r  ^r^r+1 


Case  (iii) .  a  ■  p.  and  b  ■  p  .  .  Then  .  ,  p.  .  ;  . 

\-'k  ^r^r+1  S+l^k+1  k-<-l^k+2 


S  -i  1  ® 

Vr  "-r^r+r 


in  P. 


Case  (iv) .  a  *  p 


^k^k+1 


and  b  »  o.  .  .  This  case  can  be  disoosed  of 

‘Vr 


as  in  Case  (iii) . 


Lemma  3.3  If  every  row  and  column  of  P  contains  atleast  two  positive 
elements,  then  P  contains  a  positive  loop. 

Proof.  It  is  easy  to  construct  an  infinite  sequence  of  positive  elements 
in  P  satisfying  the  hypothesis  of  Lemma  3.2. 


We  now  give  five  different  ways  of  characterizing  the  extreme 

points  of  M (?,/?-»..•  ,P  —  »q  )• 

12  tu  1  2  n 

Theorem  3.3  Let  P  be  any  matrix  in  M (o,  .p., , . . .  ,p  ;q.  ,q.. ,  • . .  ,q  ).  The 
—  ■  —  ■  ■  1  i  m  1  2  n 

following  statements  are  equivalent. 


(i)  P  is  not  an  extreme  point  of  A<(p,  ,0-, . . .  ,p  fq,  »q.,  f  •  •  •  )• 

12  m  1  2  n 

(ii)  There  is  a  positive  loop  in  P. 

(iiil  There  exists  a  submatrix  E  of  P  having  the  property  that 
every  row  and  column  of  E  has  atleast  two  positive  elements. 

(iv)  There  exists  a  square  submatrix  D  of  P  having  the  property 
that  every  row  and  column  of  D  contains  atleast  two  positive  elements. 


'j 

v\ 
•  -  •  ■  » 


(v)  There  exists  a  square  submatrix  F  of  P,  say,  of  order  k, 
having  the  property  that  the  number  of  positive  elements  in  F  is  atleast  2k.  .-.  ‘i- 


Proof .  The  equivalence  of  (i)  and  (v)  was  proved  by  Lindenstrauss  (1965, 
p.382).  The  equivalence  of  (i)  and  (ii)  was  proved  by  Klee  and  Witzgall 
(1968,  Theorem  4,  p.265).  We  prove  the  equivalence  of  these  five  statements 
directly  as  follows.  (Lemma  3.2  plays  a  crucial  role  in  our  proofs.) 


^  «  ft  %  •*, 


(i)  ««■>  (ii) .  Suppose  P  is  not  an  extreme  point  of  M(pj^,P2f •  •  • 

q  ).  For  the  remainder  of  this  proof,  we  denote  the  above  conroact 

1  2  n 

convex  set  simply  by  W.  Then  we  can  write  P  »  aB  +  (l-a)C  for  some 

distinct  B  »  (b.  .)  and  C  *  (c.  .)  in  M  and  O  <  a  <  1.  If  p  *0, 
ij  1-3  ^3 

then'  *  0  »  c^^.  Since  P  and  B  are  distinct,  there  exists  a  pair 

(i  ,3  )  such  that  o.  .  b.  .  .  Cleeurly,  p.  .  is  positive.  Assume, 

without  loss  of  generality,  that  P,  •  >  b.  .  .  There  exists  30  3i  such 

1^1  ^1^1 

that  p,  .  <  b.  .  .  Note  that  o,  .  is  positive.  There  exists  i_  it  i, 

^1^2  '■1^2  ^  ^ 

such  that  p.  .  ^  b.  .  ,  eund  3,  3o  such  that  p.  .  <  b.  and  so  on. 

^2^2  ^2^2  ^  2  ^2^3  ^2^3 

Thus  we  can  construct  a  sequence  p.  j  »  ^  ?  P.  ^  »  P.  ^  ?  .••  of 

^1^1  “1^2  ^2^2  ^2^3 

positive  elements  in  ?  with  the  property  that  any  two  consecutive  suffixes 
are  distinct.  By  Lemma  3.2,  we  can  find  a  positive  loop  in  P. 


(ii)  *=*>  (v) .  Suppose  P  contains  a  positive  loop  P.  a  ,  P^  a  ,  ...  ,  V 

^1^1  ^2^2 


^  •  Let  «  the  number  of  distinct  elements  among  i^p  i^ 


f  •  •  •  t  i 


...... ... , . ^ 

without  loss  of  generality,  that  s^  —  *2'  ^  submatrix  of 

order  s^  x  determined  by  i^^,  i^,  ...  ,  i^  -  th  rows  and  j^,  j^,  ...  ,  ' 

-th  columns.  Let  D  be  a  submatrix  of  E  obtained  by  deleting  s^  -  s^^ 

columns  from  E.  Since  every  column  of  D  contains  atleast  two  positive  ^ 

*% 

elements,  the  number  of  positive  elments  in  D  is  atleast  2s.,  and  the  order  'I'' 

1 

* 

of  the  matrix  D  is  s,  x  s, ,  ^ 

II  r 


-.'I 


(V) 
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>  (iv).  This  is  clear. 

(iv)  mmm>  (iii) .  This  is  also  clear. 

(iii)  n«>  (ii)  .  One  can  easily  extract  a  positive  loop  frota  E  and 
hence  from  P. 

(ii)  —»>  (i).  Suppose  ^  ^  '  •  •  •  '  *  positive 

1^1  2^2  k-'k 

loop  in  P.  Let  0  <  £  <  minimuni{  p.  .  ,  p.  .  ,  ...  r  p.  j  }•  Define  two 

^1^1  ^2^2  K^k 


matrices 

B  -  (b^^l 

and 

C  - 

(Ci^) 

as  follows. 

‘’l  3  ■ 

S^S 

•°S^s 

+  £ 

,  if 

s  is  odd. 

m 

-  £ 

B  if 

s  is  even, 

b  * 

ij 

"ij  ' 

if 

(i,j) 

(i.,j„)  for  any  s  »  l,2,...,k, 
s  s 

c  • 

i  j 

S''S 

"i  j 

s  s 

+  £ 

B  if 

s  is  even, 

- 

"i  j 

s  s 

-  £ 

,  if 

s  is  odd , 

c  ■ 

ij 

"ij' 

if 

(i,j)  J*  *  *  l,2,...,k. 

Note  that 

3  and 

C  are 

in 

M,  distinct  and  P  «  ‘j(B  +  C).  Consequently, 

P  is  not  an  extreme  point  of  M. 

This  completes  the  proof  of  the  theorem. 

The  above  theorem  provides  some  information  on  the  supports  of 
extreme  point 

of  the  ^bivariate  distributions  in  W(p,  »P- # . . .  rP  ;<3,  — »<?  )• 

^  1  ^  m  1  ^  n 


Corollary  3.1  Let  P  and  Q  be  two  distinct  extreme  points  of 

M (p, ,p_, . . . ,p  ;q, ).  Then  the  suooort  of  P  neither  contains 
12  m  1  2  n 

nor  is  contained  in  the  support  of  Q  properly. 


Proof.  Suppose  the  support  of  P  is  contained  properly  in  the  support  of 


Q.  Let 

be  such  that  q,  >  0  * 

*^1^1 

-°i  j  ' 
1^1 

where 

P 

*  (o  ) 

and 

Q  - 

) .  There 

exists  ^2  ^  ^1  that 

q.  • 

p.  .  . 
^1^2 

Clearly, 

is  positive.  For,  otherwise,  p  . 

^1^2 

will  be 

equal 

to 

zero.  There 

exists 

such  that  q.  .  >  o,  .  which  imolies  that 

is 

positive.  Proceeding  as  above,  we  can  construct  an  infinite  sequence 


q.  .  ,  q.  .  ;  q,  .  ,  Q;  •  ,  •••  of  positive  elements  in  Q  with  the 

^1^1  ^1^2  ^2^2  ^2^3 

property  that  any  two  consecutive  suffixes  are  distinct.  By  Lemma  3.2,  there 


exists  a  positive  loop  giving  rise  to  a  contradiction  to  the  veracity  of 


Theorem  3.3.  A  similar  argument  shows  that  the  support  of  Q  is  not 


contained  in  the  support  of  P  properly. 


Corollary  3.2  Let  P  ■  2  “  distinct  extreme  points 

of  M(p, ,p-, . - . ,p  ;q, ,q-,...,q  ).  Then  the  supocrts  of  P  and  Q  are 
12  m  1  2  n 

distinct. 


that 

Proof .  Suppose^the  supports  of  P  and  Q  are  equal.  Since  P  and  Q 

are  distinct,  there  exists  a  pair  (i.  ,j.)  such  that  P.  q,  ^  • 

1  1  *^1^1  ^1^1 


Assume,  without  loss  of  generality,  that  p  . 


which  obviously 


inplies  that  that  ?  .  is  positive.  There  exists  ^2 

p  ^  <5.  •  Since  the  supports  are  equal,  p  .  is  positive.  There 

V2  "-1-2  1^2 

exists  i  i.  such  that  o.  .  ^  q.  .  ,  and  so  on.  Thus  we  can  construct 
2  1  "^2^2  '^2^2 

an  infinite  sequence  P, .  ^  »  P.  ^  ?  P,  ^  »  P,  4  f  positive 

1^2  2^2  2-'3 

elements  in  P  with  the  property  that  any  two  consecutive  suffixes  are 
distinct.  By  Lemma  3.2,  there  exists  a  positive  loop  in  P  negating  the 
(validity  of  Theorem  3.3.  This  contradiction  proves  the  result. 


Corollary  3.3  Ihe  number  of  extreme  ooints  of  M(p, ,p,,. . . ,p  ;q, ,q-,. . . ,q  ) 
■■  •  X  ^  m  X  ^  n 

is  finite. 

Proof.  The  set  of  all  subsets  of  { (i, j)  ;  r  *  1  to  m,  j  =  1  to  n}  is 
finite. 


corollary  3.4  Let  S  » 


where  s^^  ■  ^ 

Mtpj^  ,?2 » •  •  • 'Pjjj'^2^  *^2 ' "  *  * ^ 


q^  ■  1  for  some  j. 


be  the  matrix  in 
and  j .  Then  S 
:  and  only  if  p^ 


M  (p^  ,©2 » •  •  •  '^2 '  * 

is  an  extreme  point  of 

*  1  for  some  i  or 


Proof .  If  either  p^  ■  1  for  some  i  or  q^  ■  1  for  some  j,  then 

W(p- ,0  jq,»q-» — )  contains  only  one  matrix  S  and  consequently, 
12  *  m  1  2  n 

it  is  an  extreme  ooint  of  M(p,  ,0-,. . .  ,p  ;q.  ,q., , . . .  ,q  ).  On  the  other  hand, 

12  m  1  2  n 

if  p^  1  for  every  i  and  q^  1  for  every  j,  then  we  can  find 

Py*  P^»  positive  with  u  v  and  r  t.  Then 


is  a  positive  loop  in  S.  Hence  S  is  not  an  extreme  point. 


The  matrix  S  defined  above  is  the  product  measure  of  the  distribn-. 

p  ,p  and  q, .q, . q  -  the  above  corollary  is  a  discrete  analogue 

*^12  o  ^  ^  “ 

of  Theorem  3  of  KeB5>  (1968,  p.l356). 

Examples. 

1.  m  -  2,  n  -  2;  ■  Pj  *  's  *  *  ^2' 

Extreme  points  of  *te 


\  o1 

b  ‘s' 

and 

P  ’ll 

js  q 

2.  m  -  2,  n  -  2;  p^  -  1/4.  P2  “  3/4;  q^  -  2/3  and  q^  -  1/3. 
Extreme  points  of  M (1/4, 3/4;  2/3, 1/3)  are 


”0 

3/12 

3/12 

0  ' 

_8/12 

1/12 

and 

_5/12 

4/12_ 

3.  m  -  2,  n  •  3;  “  ®2  *  “  ^2  "  **3  * 

Extreme  points  of  M (1/2, 1/2;  1/3, 1/3, 1/3)  are 


"2/6  1/6 

0  ' 

"2/6 

0 

1/6" 

0  1/6 

2/6_ 

,  [o 

2/6 

1-/6. 
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"  0 

1/6 

2/6" 

and 

1/6 

0 

2/6“ 

_2/6 

1/6 

0 

1/6 

2/6 

0 

Extreme  points  of  M _ (P, *P^» — »P_J  » • • • fO. ) 

— — —  —  pop  X  ^  m  1  ^  n 

We  discuss  a  method  of  enumerating  the  extreme  points  of  the 
above  conmact  convex  set  with  the  help  of  some  examples.  From  this  discussion, 
a  general  strategy  can  be  devised  to  enumerate  the  extreme  points  in  the 
general  case. 

1 .  The  2x2  case. 

One  can  show  that  any  bivariate  distribution  in  MpQjj(Pj^*?2' 
given  in  the  form  of  a  matrix  P  «  must  satisfy  the  inequality 

i  P22  i  '’2''  '52- 

where  a  A  b  denotes  the  minimum  of  the  two  numbers  a  and  b.  Conversely, 
if  a  number  p^^  satisfies  the  above  inequalities,  one  can  find  a  matrix 
in  MpgQ(P2^»P2'  whose  (2,2)-th  element  is  p^^*  More  precisely, 

this  matrix  is  given  by 


Pi-<?2‘^®22  '^2'’°22 


P2-P22 


Consequently,  the  extreme  points  of  »P2»’  obtained  by 

setting  •  P2'^2  ^22  *  ^2  ^  ^2  rest  of 

the  entries  of  the  matrix  using  marginality  conditions.  There  are  only  two 


y.v..-.y. 


■  **  •**  «'•  ,*•  *■ 


extreme  points  of  '  pQD^^l'°2'  these  are  given  by 


Piqi  ?^q2 


•Vl  -°2‘^2 


Piqi  Piq2j 


P2qi  P2<J2 


q2-P2 


*°2-^2  ‘’2 


if  P2  ^  I2 


if  P2  A  ‘I2 


Every  member  of  '°2'  *  convex  combination  of  these 

two  extreme  points.  The  following  are  two  concrete  examples. 


Example  1.  ^1  “  *  *^1  *  ^2' 


The  extreme  points  are 

1/2  0” 

and 

_0  1/2 


1/4  1/4“ 

1/4  1/4 


Example  2.  p^^  ■  1/4,  p^  ■  3/4  and  ■  2/3,  »  1/3. 


The  extreme  points  are 


0 

and 

"2/12 

1/12 

Ls/12 

4/12 

_6/12 

3/12_ 

2.  Ihe  2x3  case. 


In  this  case,  the  determination  of  the  extreme  points  of 


fP-,'*  'll '‘I-j)  can  be  achieved  graphically.  If  the  matrix  (p..) 

PQu  1  2  1  2  J  i  3 

belongs  to  ^pQo^P2.'^2''  *51  *'^2 '^3^'  check  that  the  following 

inequalities  must  be  satisfied. 

P2'53  i  ^'23  i  •‘'2  ^“*3  “■* 

(P2'2j  ^  V  Pj3  i  Pj2  *  Pj3  i  =2  A  (qj  *  PJ3I. 

(a  W  b  means  the  maximum  of  the  two  numbers  a  and  b.)  Conversely,  if 
P22  and  P22  are  two  nvimbers  satisfying  the  above  inequalities,  then 
the  matrix 

qi-P2^22*P23  ^^22  *^3-^23 

^2’°22"°23  “22  °23 

belongs  to  ^pqd^^1'^2^  ^1'‘*2'^3^'  above  two  inequalities  determine 
a  simplex  in  the  two  dimensional  P22  “  P23  plane.  The  extreme  points  of 
this  simplex  give  the  extreme  points  of  '^2'  °l''^2'^3^’  ^ 

illustration,  let  °i  "  ^2  "  *^1  *  *^2  *  ^3  *  determining 

inequalities  are 

1/6  P22  £  1/3  and 

(1/3)  V  P23  “  1/3  <  P22  +  ?23  1  1/2  A  (1/3  +  P23)  -  1/2. 

These  inequalities  determine  the  following  simplex  in  the  p^^  ~  P23  plane. 

There  are  four  extreme  points  of  the  set  M  Cl/2,1/2;  1/3, 1/3, 1/3)  which 

?QD 

are  determined  by  the  four  extreme  points  of  the  simplex.  Note  that  the  set 
M(l/2,l/2;  1/3, 1/3, 1/3)  has  six  extreme  points. 


The  four  extreme  points  corresponding  to  P^,  ^3*  **4  respectively  are 


1/6 

1/6 

1/^ 

"1/6 

1/3 

o' 

1/6 

1/6 

1/6^  , 

_l/6 

0 

l/3_ 

1/3 

1/6 

0 1  ri/3 

0 

1/6 

_0 

1/6 

1/3J  ,  \_0 

1/3 

l/6_ 

Every  member  of  M  «(l/2,l/2;  1/3, 1/3, 1/3)  is  a  convex  combination  of  these 
PQD 

four  extreme  points. 

3.  The  2x4  case. 

One  can  check  that  if  a  matrix  belongs  to  ^pqd^P]^'P2' 


-='2'’4  i  -°24  i 


'^2''3  ^2'‘4'  «  ^24  -  ^23  "24  -  *'2  ”  ”3  "24'  ’ 

(P2'J2  *  P2'=3  ■"  ‘’2'’4>  V  <P23  *  ^’24’  i  *’22  *  ''23  •"  '’24  i 
'2  ‘’2  *  '23  ■'■''24>- 

conversely,  if  ^22’  ^23  °24  three  numbers  satisfying  the  above 

inequalities,  then  they  determine  a  member  of  ^pqd^P2.'°2'  '^1'^2'^3^ 
an  obvious  way.  For  details,  see  Subramanyaua  and  Bhaskara  Rao  (1985) .  We 
illustrate  the  determination  of  the  extreme  points  of  this  set  with  the 
help  of  a  concrete  example. 

“  ?2  *  •  <12  *  %  “  *^4  *  determining 

ineqtialities  are 


1/8  <  Pj^  <  2/8, 

2/8  1  P23  P24  1  2/8  +  p^^  and 

3/8  V  (P23  +  P24)  <  P22  +  P23  +  P24  -  ^  ?23  ^24’- 

The  first  step  in  the  determination  of  all  extreme  points  consists  of 
eliminating  the  symbols  V  and  A  from  the  above  set  of  inequalities 
by  splitting,  if  necessary,  some  of  the  inequalities.  For  example,  the 
inequality  2/8  £  p^^  +  Pj^  ^  2/8  +  p^^  is  equivalent  to  the  two 

inequalities  2/8  <  p_-  +  p_ .  <  3/8  and  3/8  <  p.,  +  p..  <  2/8  +  p 
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The  splitting  certainly  helps  to  get  rid  of  the  syinbols  V  and  A  .  The 
above  set  of  determining  inequalities  is  equivalent  to  the  following  two 
sets  of  ineoualities. 


1/8 

<  *0 
-  -24 

<  2/8 

2/8 

^  -“23 

^P24  1 

3/8 

3/8 

-  -22 

^  P23  +  P 

24  i 

1/8 

1  P24 

<  2/8 

3/8 

-  ^23 

+  < 
^24  - 

2/8  .  p^. 

o  +o  <o  +D  +0  <  4/8 

*23  ^24  -  ^22  ^23  ^24  - 

Now,  the  second  step  consists  of  the  following  mancevres.  In  each  set 

of  inequalities,  set  the  central  expression  equal  to  the  expression  either 

on  the  right  or  on  the  left.  This  will  result  in  three  equations  in  three 

unknowns.  These  linear  equations  are  very  easy  to  solve.  From  the  solution 

so  obtained,  one  can  build  a  matrix  in  Mp^^(l/2,l/2;  1/4 ,1/4 ,1/4 ,1/4) 

using  the  marginality  restrictions.  The  set  of  all  extreme  points  of 

Mp^jj(l/2,l/2;  1/4, 1/4, 1/4, 1/4)  is  contained  in  the  collection  of  all 

matrices  so  obtained  above.  After  eliminating  the  duplicates,  the  following 
the 

is  collection  of  all  extreme  points  of  M^^^ll/2,1/2;  1/4, 1/4, 1/4, 1/4) . 


‘1/8 

1/8 

1/8 

1/87  f2/8 
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1/8 

1/8 

1/^  ,  |_  0 

2/8 

1/8 

l/8_ 

illeations . to  Continggncy  Tables 


In  this  section,  we  are  Interested  in  the  following  problem.  Let 
X  and  Y  be  random  variables  with  known  marginal  distributions  X  and  v 
respectively  but  with  unknown  joint  distribution  function.  We  want  to 
test  the  hypothesis  that  X  and  Y  are  independent  against  the  alternative 
that  X  and  Y  are  strictly  positive  quadrant  dependent.  For  sin^licity. 


A  problem  of  choice  arises  to  test  the  above  hypothesis  of  independence 


against  the  above  specific  alternative.  One  way  to  resolve  the  dilemma  is 


is  M  (X,\)).  Conputation  of  the  newer  of  a  given  test  at  every  point 
PQD 

in  the  domain  is  not  feasible  practically.  The  following  theorem  asserts 
that  it  is  enough  to  compute  the  power  at  all  extreme  points  of 

Theorem  4.1  Let  the  extreme  points  of  •••  » 

.  Let  T  be  a  test  proposed  to  test  the  hypothesis  of  independence 

and  S«(  )  its  newer  function.  Let  U  be  an  arbitrary  distribution 
T 

in  •  (Then  we  can  write 

PQD 


...  dj^y 


for  some  ^^2'  ***  ' 


>  O  with  ®2  “*  \ 


B^(y)  »  +  \ 


For  a  discussion  of  this  theorem,  see  Bhaskara  Rao»  Krlshnaiah  and 
(1985). 

Subramanyam  >1  In  view  of  this  theorem,  if  we  wish  to  compare  the  performance 
of  two  tests,  we  merely  compute  the  powers  of  these  tests  at  the  extreme 
point  distributions  and  then  coispare  these  powers  point  by  point.  To 
illustrate  the  mechanism  of  this  theorem,  we  consider  the  case  m  >  2  *  n. 


Asstane,  for  sinplicity,  that  p^  £  <12’  extreme  points 


There  are  two  popular  measures  of  association  between  X  and  Y 


defined  by 


and 


Gamma  Ratio 


Y 


^11^22  ~  ^12^21 
^11-22  *  °12^21 


(see  Goodman  and  Kniskal) 


P 

Spearman's  rho  *  p  *  — 

If  X  and  Y  are  independent,  then 
show  that  Y  >  0  if  and  only  if  X 


11^22  ~  ^12^21 

Y  •  O  *  p.  Further,  one  can 
and  Y  ar^  strictly  positive 


quadrant  dependent  and  that  p  >  0  if  and  only  if  X  and  Y  are  strictly 


positive  quadrant  dependent. 


Test  based  on  the  ganma  ratio  ;  T^^ 


Reject  the  null  hypothesis  if  and  only  if 

Y  •  ~  "l2”21  ^ 

"ll”22  “l2"21 

Test  based  on  the  rho  criterion  :  T^ 

Reject  the  null  hypothesis  if  and  only  if 


Let  us  calculate  the  power  of  the  tests  and  at  the 

extreme  points  of  ^pqd^^1'^2'  °1'‘*2^‘ 

(P, )  ■  S_  (P, )  »  Sire  of  the  tests  T,  and  T, . 

1  T2  1  12 

A 

Under  P^,  «  O  almost  surely  and  consequently,  y  “  1 

with  probeibility  unity.  Therefore, 

S  (p  )  ■  Pr (Rejecting  the  null  hypothesis/P.) 

2 

A 

«  Pr(y  >  a/  P^)  -  1 

It  is  now  obvious  that  8_  (P_)  <  6_  (P_) •  Hence  the  power 

T2  2  “  2 

function  of  dominates  the  power  function  of  Tj.  As  a  matter  of  fnctr 

the  power  function  of  7^  will  always  dominate  the  power  function  of  any 
test  proposed  for  testing  the  hypothesis  of  independence.  In  a  nutshell, 
what  this  means  is  that  the  test  based  on  Gamna  Ratio 

is  the  uniformly  most  powerful  test  for  testing  the  hypothesis  of 
independence  against  the  alternative  of  strict  positive  quadrant  dependence 
in  the  context  of  2x2  contingency  tables. 

Conraarisons  of  the  power  functions  of  various  tests  have 

been  carried  out  more  eladwrately  in  Bhaskara  Rao  ,  Krisnnaiah  and 
Subromanyam  ( 1985 '  • 
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